DEUTSCHES ELEKTRONEN-SYNCHROTRON 

DESY 02-048 
May 2002 



Transverse self-fields within an electron bunch moving 

in an arc of a circle 



Gianluca Geloni, Jan Botman, Jom Luiten and Marnix van der Wiel 

Department of Applied Physics, Technische Universiteit Eindhoven, 
P.O. Box 513, 5600MB Eindhoven, The Netherlands 



Martin Dohlus, Evgeni Saldin and Evgeni Schneidmiller 

Deutsches Elektronen-Synchrotron DESY, 
Notkestrasse 85, 22607 Hamburg, Germany 



Mikhail Yurkov 

Particle Physics Laboratory (LSVE), Joint Institute for Nuclear Research, 
141980 Dubna, Moscow Region, Russia 



Transverse self-fields within an electron bunch moving in an arc of a circle 



Gianluca Gelonijj Jan Botman, Jom Luiten, and Marnix van der Wiel 
Department of Applied Physics, Technische Universiteit Eindhoven, 
P.O. Box 513, 5600MB Eindhoven, The Netherlands 

Martin Dohlus, Evgeni Saldin, and Evgeni Schneidmillcr 

Deutsches Elektronen-Synchrotron DESY, 
Notkestrasse 85, 22607 Hamburg, Germany 

Mikhail Yurkov 

Particle Physics Laboratory (LSVE), 
Joint Institute for Nuclear Research, 
141980 Dubna, Moscow Region, Russia 

As a consequence of motions driven by external forces, self-fields (which are different from the 
static case) originate within an electron bunch. In the case of magnetic external forces acting 
on an ultrarelativistic beam, the longitudinal self-interactions are responsible for CSR (Coherent 
Synchrotron Radiation)-related phenomena, which have been studied extensively. On the other 
hand, transverse self-interactions are present too. At the time being, existing theoretical analysis 
of transverse self-forces deal with the case of a bunch moving along a circular orbit only, without 
considering the situation of a bending magnet with a finite length. In this paper we propose an 
electrodynamical analysis of transverse self-fields which originate, at the position of a test particle, 
from an ultrarelativistic electron bunch moving in an arc of a circle. The problem will be first 
addressed within a two-particle system. We then extend our consideration to a line bunch with a 
stepped density distribution, a situation which can be easily generalized to the case of an arbitrary 
density distribution. Our approach turns out to be also useful in order to get a better insight in the 
physics involved in the case of simple circular motion and in order to address the well known issue 
of the partial compensation of transverse self-force. 

PACS numbers: 



I. INTRODUCTION 

When an electron bunch undergoes a motion under the 
influence of external forces, the particles within the bunch 
become sources of self-fields which are different from the 
static case and which obey the usual Lienard-Wiechert 
expressions. 

These self- fields feed back on the particle dynamics, 
which often makes the description of the evolution of the 
system a problem without easy solution. 

Many electrodynamical systems fit in the latter de- 
scription. In all these systems one can recognize two 
separate aspects in the evolution issue: a dynamical one, 
which is governed by the equation of motion, and an 
electrodynamical one, which is taken care of by Maxwell 
equations. The solution to the evolution problem is ob- 
tained when one is able to solve simultaneously both 
equations: this is, for example, what self-consistent com- 
puter codes strive for. An example of such self-consistent, 
numerical solutions is given (see Q|) by the program 
TRAFIC 4 (which will be employed in this paper to 
provide cross-checks with analytical results). Although 
TRAFIC 4 and other codes are able to provide quite re- 
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liable results for a number of evolution problems, a sep- 
arate study of their dynamical and electrodynamical as- 
pects, which can be achieved by means of theoretical 
analysis, is necessary in order to get a good understand- 
ing of the physics involved, and to provide benchmark 
and cross-checks for simulations (and, at an earlier stage, 
to build correct simulations too). 

Moreover, it often happens that this kind of theoreti- 
cal analysis can be directly applied in order to get good 
practical solutions for the problem at least in a narrow 
region of the parameters which specify the system setup. 
The results can, then, be used for quick estimations of 
the magnitude of the effects under investigation. 

Let us restrict our attention to the case of self-fields 
within an ultrarelativistic electron bunch driven by ex- 
ternal forces of magnetic nature. 

The self-interaction in the longitudinal direction (par- 
allel, at any time, to the velocity vector by definition) 
is responsible for the energy exchange between the sys- 
tem and the acceleration field and for all CSR (Coherent 
Synchrotron Radiation)-related phenomena, which have 
been studied extensively elsewhere (see, among the oth- 
ers, §... §). ' 

These investigations are important in view of the need 
for very high-peak current, low emittance beams for self- 
amplified spontaneous emission (SASE)-free-electron- 
lasers operating in the x-ray region (see, among the oth- 
ers, P|). Similar beams are also being considered for 
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production of femtosecond radiation pulses by simpler 
schemes based on Cherenkov and Transition Radiation 
JTo| ]: production and utilization of such kind of beams 
may prove difficult due to self-field collective effects. In- 
deed, these effects may spoil the required high brightness 
of the electron beam, which is a matter of major concern 
among people involved this kind of physics research. 

The study of self-forces in the transverse direction is 
important for the same reasons. They were first ad- 
dressed, in the case of a circular motion, and from an 
electrodynamical viewpoint, in |fj"l[ . Further analysis ((g), 
p2[... p^|) consider, again, the case of circular motion 
both from an electrodynamical and a dynamical view- 
point and in the approximation of a rigid bunch: at the 
time being, existing theoretical analysis of transverse self- 
forces deal with the case of circular orbit only, without 
considering transient collective phenomena. 

In this paper we propose a fully electrodynamical anal- 
ysis of transverse self-fields originating, at the position of 
a test particle, from an ultrarelativistic electron bunch 
moving in an arc of a circle, thus treating for the first 
time, besides the basic situation of circular motion, also 
the case of transient between a straight line and a hard- 
edge magnet (and, vice versa, from a bend to a straight 
line) . 

Consistently with the choice to analyze the electrody- 
namical aspect of the problem only, a zero energy-spread 
will be understood when considering the evolution of an 
electron bunch. As underlined before, although the re- 
sults obtained can be directly applied, from a practical 
viewpoint, only in the case in which the zero energy- 
spread hypothesis is verified a posteriori, the outcomes 
of this paper are important to get a good insight if the 
physics of the problem, and to provide benchmark and 
cross-checks for simulations. 

Firstly, a two-particle model is adopted in order to 
study the transverse force produced by a single parti- 
cle, and then, by summing up all the contributions from 
different electrons, the case of a line-bunch model charac- 
terized by a rectangular density distribution is analyzed: 
this can, in fact, be easily generalized to the case of an 
arbitrary density distribution. 

Besides providing results in the case of a finite hard- 
edge bending magnet, our approach turns also useful in 
order to get a better insight of the physics involved in 
the case of simple circular motion. 

Moreover, our results give us a better understanding 
(although one has to take into account, here, dynamical 
aspects as well as electrodynamical ones) of the partial 
compensation between transverse self-force and gradient 
of the potential energy deviation from the nominal value. 

We will discuss the limits of applicability of our model, 
with respect to the transverse beam size, h, in Section 
III). Within such applicability limits, we find results 
which are independent of the bunch transverse dimen- 
sion. 

The paper is organized as follows. The transverse in- 
teraction between two electrons moving on a circle, to- 



gether with a simple dyna mic al interpretation is treated 
in Section H. In Section III, we deal with a stepped- 



profile electron bunch interacting with a test particle 
again on a circle, and discuss also the applicability re- 
gion of the line model. Transient behavior (from straight 
to circular path and vice versa) for the transverse self- 
forces between two particles are then studied in Section 



IV. Results for the transient of a stepped-profile bunch 
are given in Section [v], where a treatment for the case of 
a more generic bunch density is also proposed. A regu- 
larization technique for cancelling the singularity in the 
expression for the transverse force (which always arises 
in the limit of a near-zero distance between test particle 
and sources) is then applied. In Section |Vl] we will deal 
with the well-known issue of the partial com pensation of 
transverse self- force. Finally, in Section VII, we come to 



a summary of the obtained results and to conclusions. 



II. TRANSVERSE INTERACTION BETWEEN 
TWO ELECTRONS MOVING IN A CIRCLE 

Let us begin our study considering the steady case of 
two electrons moving on a circle of radius R. The electro- 
magnetic force which one of the two particles (designated 
with "T", i.e. the test particle) feels, due to the inter- 
action with the other one (designated with "S", i.e. the 
source particle), is given by 



F(r T ,t) = eE(r T ,t) + ec(3 T x B(r T ,t), 



(1) 



where rx is the position of the test particle, e is the 
electron charge with its own (negative) sign, /3x is the 
velocity of the test particle normalized to the speed of 
light, c, while E(rx,t) and B(rT,t) are, respectively, 
the electric and the magnetic field generated at a given 
time t by the source particle S, at the position of the test 
particle T, namely 



E{r Tl t) 



(n - /3 S ) x fa 



and 



R ST (1 - fi • AO 



B(r T ,t) = -fix E(r T ,t) . 

c 



(2) 



(3) 



Here (3s and 0s are, respectively, the dimensionless ve- 
locity and its time derivative at the retarded time t' , 
Rst is the distance between the retarded position of the 
source particle and the present position of the test elec- 
tron, n is a unit vector along the line connecting those 
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two points and 7s is the usual Lorentz factor referred to 
the source particle at the retarded time t'. 

The transverse direction (on the orbital plane) is, by 
definition, orthogonal to (3t. The transverse component 
of Eq. (0) can be written as the sum of contributions 
from the velocity ("C", Coulomb) and the acceleration 
("R", Radiation) fields, namely 



F±(rr, t) = Fj_ c (r T , t) + Fj_ R (r T , t), 



(4) 



where 



and 



(T,f) 



n x (1 - /3 S • fa) ~ Pxs (1 - & • Pt) 



7 S 2 R 2 ST (1 - n • PsY 



(5) 



_LR 



(rr,t) 



4tt£o 



n x (n • p s ) (l-Ps- Pt) - P±s (n • p s ) (1 - n • /3 T ) /3xt (1 - n • /3t) + n± (#r ■ /9s 



ii ST (1 - fi • p S ) 



R ST (l-n-Ps) 




FIG. 1: Geometry for the two-particle system in the steady 
state situation, with the test particle ahead of the source. 
Here T is the present position of the test particle, S is the 
present position of the source, while S' indicates the retarded 
position of the source. 



F ±c 



1 + /3 2 -2/?cos((/>/2) 



47T£ 4i?V sin(0/2)(l -/3cos(</>/2)) 3 



3 2 



47T£ 2i? 2 



sin(0/2)(l-/3cos(</>/2)) 



(1 -p 2 ) sin (f>/2 



(l-/?cos(0/2)) d 



(7) 



(8) 



If we now assume a small retarded angle < 1), which 
is justified in the ultrarelativistic case, we can expand 
Eq. (^) and Eq. (||) to the second non-vanishing order 
in <fi thus obtaining 



F ±c 



e 2 7 3 

47T£ i? 2 



(9) 



and 



Let us first consider the case in which the test particle is 
in front of the source. In this case, referring to Fig. |l|, we 
can define with As the curvilinear distance between the 
present position of the test and of the source particle; 4> 
will indicate, instead, the angular distance between the 
retarded position of the source and the present position of 
the test electron, and it will be designated as the retarded 
angle. 

In the following we will assume (3s — Pt — ft- This hy- 
pothesis will naturally lead, further on, to the assumption 
of a zero energy-spread when one considers the evolution 
of an electron bunch, and, as already discussed in Sec. §, 
it is consistent with our choice to analyze only the elec- 
trodynamical aspect of the problem. Upon this we can 
write Eq. (^|) and Eq. (^) in the following way: 



Fx 



e 2 7 3 



where we define $c an d as 



*r(0) , 



$cW>) 



(1 + 2 /4)3 



and 



0(l + 2 /4) 3 



(10) 



(11) 



(12) 



Here and above <fi = 7^. This normalization choice, al- 
ready treated in |3| , is quite natural, 1 being the syn- 
chrotron radiation formation-angle at the critical wave- 
length. In the derivation of Eq. (|^) and Eq. (|l0|) (and 



in the following, too) we understood <j> 3> 1/7, which 
is again justified by the ultrarelativistic approximation. 
Moreover, it is important to realize that the assumption 
(ft *C 1 is by no means a restrictive one, since it keeps 
open the possibility of comparing (ft with the formation 
angle 1/7 (note that a deflection angle smaller or bigger 
than I/7 is characteristic of the cases, respectively, of 
undulator or synchrotron radiation). 

The following expression can be then trivially derived, 
which is valid for the total transverse force felt by the 
test particle 



<E>(As)As 



where <f> is defined by 



e 2 7 3 
Ane R 2 



(13) 



H<f>) 



(1 + </, 2 /4) 3 



(14) 



Note that Eq. (|1J) is completely independent of the 
parameters of the system: it is then straightforward to 
study the asymptotic behaviors of $. In order to do so, 
just remember that the retardation condition linking As 
and (ft is given by (see 0) 



As = R<j)-2l3Rsm^ , 



(15) 



or its approximated form 



As = (1 - f3)R(f> - 



R(ft 3 
~24~ 



(16) 



It is now evident that <&(As) — > 1/(3 As) when (ft ^> 1 
and $(As) — » l/(As) when (ft <C 1, having introduced 
the normalized quantity As = (7 3 /i?)As. This normal- 
ization choice is linked to the fact that the critical syn- 
chrotron radiation wavelength, -R/7 3 , is also the minimal 
characteristic distance of our system: two particles nearer 
than such a distance can be considered as a single one ra- 
diating, up to the critical frequency, with charge 2e (see 
§)■ ' 

The asymptotic behavior above suggests to study the 
function 4>(As)As. We plotted such a function in Fig. || 
(together with the radiative contribution alone) for val- 
ues of As running from to 5. 

It is interesting to underline the fact that, as one can 
see from Fig. |^, the transverse force is always centrifugal, 
for any distance between the two particles. This fact can 
be explained by means of a simple relativistic argument 
which holds, qualitatively, for all particle distances: in 
order to build the two-particle system, i.e. to bring them 
together, one needs to work against the electromagnetic 
field. Then, the total mass of the system accounts for 
this interaction energy too, and is therefore bigger than 
the simple sum of the particles masses. Hence, also the 
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FIG. 2: Plot of $(As)As (solid line) and comparison with 
the asymptotic values, 1 and 1/3. The dashed line shows the 
radiative contribution <E>r(As)As alone. 



equilibrium orbit radius must be bigger than R, and a 
centrifugal self-force is to be expected. 

It is also worthwhile to note that <E>r has the same 
asymptotic behavior as and that <I>c gives impor- 
tant contributions to $ only for values between the two 
asymptotes (see Fig. |J). This may find intuitive expla- 
nation in the following reasoning. As is well known (see, 
for example, p5[), the velocity field of an ultrarelativistic 
electron is radial, with respect to the "virtual" position 
which the particle would assume if it moved with con- 
stant velocity starting from the retarded point, but the 
line forces are not isotropically distributed, and resemble 
more and more the plane wave configuration as j3 — > c. 
Therefore, the test particle is influenced by the velocity 
field of the source particle only when such a field " shines 
right on the test electron" (quoted from Q), which does 
not happen for asymptotic values of (ft. 

Let us now analyze, in the framework of our line model, 
the case in which the source particle is in front of the test 
^article. The geometry is qualitatively sketched in Fig. 

The difference with respect to the situation in which 
the test electron is in front of the source is that the test 
electron "runs against" the electromagnetic signal emit- 
ted by the source, while in the other case it just "runs 
away" from it. Therefore the relative velocity between 
the signal and the test electron is equal to (1 + (3)c, in- 
stead of (1 — (3)c in the other situation. Hence the retar- 
dation condition reads 



As = Rcft + 2(3Rsm- 
or, solved for (ft in its approximated form, 

As 



R(l + 0) 



(17) 



(18) 
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FIG. 3: Geometry for the two-particle system in the steady 
state situation, with the source particle ahead of the test one. 
Here T is the present position of the test particle, S is the 
present position of the source, while S' indicates the retarded 
position of the source. 



In this situation, (3s is almost parallel to (and equal to) 
/3t and antiparallel to h: it turns out that the only im- 
portant contribution to the transverse force is given by 
the second term on the right side of Eq. (||), and it is 
easy to check that 



and y, individuated by the tensor indexes 1 and 2), thus 
giving a vanishing integral contribution. All the quanti- 
ties between round brackets are calculated in the system 
rest frame. 

On the other hand, the total energy of the system turns 
out to be (see, again, [la]): 



(21) 



U syst = 7 2mc^ + T 33 dV 



where, again, the quantities between round brackets are 
calculated in the system rest frame. 

It is now easy to calculate U and T 33 in the rest frame. 
The electromagnetic interaction energy is given by the 
work done against the field to bring the two particles 
together from a situation in which they are separated by 
an infinite distance: 



U : 



47T£o7As 



(22) 



In order to calculate T33, we can consider the electrostatic 
stress tensor alone because in the rest frame, at short 
distance 7 As <C R/j 2 , the radiative field contributions 
are unimportant and, on the other hand, the external 
magnetic field has zero component along the direction of 
motion. Now, the electrostatic stress tensor is given by 



F± 



47T£ i?As 



(19) 



It may be worthwhile to underline that the force in Eq. 
(@) is, evidently, of the same magnitude as the force in 
Eq. H). 

As a note to the entire steady-state case discussion, 
let us provide an interpretation for the electromagnetic 
transverse forces from the viewpoint of relativistic dy- 
namics in the limit As <C R/j 3 - In such a situation the 
two electrons are near enough to consider them travelling 
with the same velocity vector. One can, then, analyze the 
situation in the rest frame of the two-particle system. By 
means of a Lorentz transformation, one can find that the 
total momentum of the system in the laboratory frame 
and in the direction of motion (let us call it z-direction) , 
is given by (see [[15)): 



P m = 7 2m + U/c 2 



T 33 dV (3c 



(20) 



where the integration (in Eq. ( |20| ) and below) is per- 
formed over all space. Here m is the electron rest mass, 
U is the electromagnetic interaction energy of the system 
and T33 is the component of the electromagnetic stress 
tensor corresponding to the flux of the z-component of 
the total momentum along the z direction. It is worth- 
while to remark that we do not have to include X13 and 
T 23 (corresponding anyway to momentum flux in the z- 
direction) in the calculation of P z because they are anti- 
symmetric in z (as well as in the other two directions, x 



T l3 = eoiEiEj - S l3 E 2 /2) , (23) 

where i, j — 1... 3. For us, the only interesting compo- 
nent is 



T33 - 

It can be proven that 



-E 2 Z - E 2 /2 



T33 dV = -U . 



(24) 



(25) 



Now, the equations for the momentum and the energy 
of the system in the laboratory frame read: 



and 



P z = 7 (2m + 2U/c 2 ) (3c 



U syat =7 [2mc 2 + U{1 + (3 2 )} 



(26) 



(27) 



From the transverse component of the equation of motion 
for the system one gets 



F ±syst ~ 2eB(3c - 2 



AneoRAs 



(28) 



which justifies Eq. ( [19| ) and the asymptotic behavior for 
short distance between the particles of Eq. (fi~3|) . 
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The discussion above underlines the fact that, in or- 
der to get a correct dynamical interpretation of the elec- 
trodynamical transverse forces on the two particles, one 
needs to account for the self-interaction energy and mo- 
mentum flux of the system. Nevertheless, by doing this, 
one can easily find that the energy of our system (divided 
by c) and its momentum do not constitute a four-vector 
anymore, as it is seen directly from Eq. (^0|) and Eq. 
( pl| ) . At the beginning of the 20th century a similar prob- 
lem was encountered by people trying to build a classical 
model of the electron: in that case, the solution was the 
introduction of the Poincare stress tensor (see, among 
others, jl5| or [^6|) , which granted stability to the sys- 
tem and recovered covariance for the energy-momentum 
pair. As said in [^5): "It is not unreasonable then to in- 
clude Poincare stresses in our classical models of charged 
particles, or at least to remember that care must be taken 
in discussion of purely electromagnetic aspects of such 
models". In our case we deal with a completely elec- 
trodynamical system, so that there is no straightforward 
way to introduce the analogue of Poincare stresses: nev- 
ertheless one must take care and remember that, within 
our accuracy, energy and momentum are no longer form- 
ing a four-vector. It might be worthwhile to underline 
the fact that, in modern physics of particle acceleration, 
problems which have been of pure academic importance 
for about one hundred years are now becoming of prac- 
tical importance. 



III. TRANSVERSE INTERACTION BETWEEN 
AN ELECTRON AND A BUNCH MOVING IN A 
CIRCLE 

In this Section we discuss the transverse force felt by 
an electron due to the interaction with a line bunch (with 
rectangular density distribution) moving in a circle. The 
geometry is described in Fig. EL and will be considered 
rigid, i.e. fixed during the entire evolution of the system. 
This fact means that we will have zero energy-spread in- 
side our bunch and this is in agreement with our model 
choice in Section || and our program of a pure electrody- 
namical analysis, without dynamical aspects. 

Before beginning actual calculations, we provide here 
a discussion about the applicability region of our line 
model. In the case of contributions from particles be- 
hind the test particle, the region of applicability of the 
line model follows straightforwardly from the retarda- 
tion condition: one can easily check that the inclusion 
of a transverse dimension of the bunch, already desig- 
nated with h in Section B, adds to Eq. ( |l6| ) the term 
of magnitude h 2 /(R(f>). Such a term is negligible with 
respect to the others in the retardation condition, when- 
ever h -C -R/7 2 , which specifies the region of applicability 
of our model as regards the transverse bunch size. 

The situation becomes more complicated when one 
considers contributions from electrons in front of the test 
particle. In fact, in the case As < (source particle in 



front) and |As| < h, we have a situation in which the 
test electron overtakes the source before it is reached by 
the electromagnetic signal: this means that, even if the 
test particle is behind the source, the present position of 
the test particle is, anyway, in front of the retarded posi- 
tion of the source. Then, it is possible to show that the 
line model constitutes a valid description of the situation 
only when |As| 3> h. The cases that do not verify such 
a condition are left for future study. 



FIG. 4: Schematic of the test electron T and the stepped- 
profile bunch; As m i n indicates the distance between the elec- 
tron and the head of the bunch; As max is the sum of the 
bunch length and As m i„ 

Now, in order to actually evaluate the transverse elec- 
tromagnetic force which the bunch exercises on the test 
particle one has to sum up the contributions from all the 
retarded sources. Since the two-particle interaction has 
been calculated in Section [n] as a function of the retarded 
angle 4>, it is convenient to switch the integration variable 
from As to <b. 



= R 1-/3 cos- 



(29) 



being the Jacobian of the transformation. Therefore, 
from Eq. (g), Eq. (§) and Eq. @ one has 



F 



_LC 



e 2 A 1 
AneoR 27 2 



F± c (<t>) i?(l-/3cos(0/2)) 

(3 



1 -/3cos(0/2) 



- ln(tan(0/4)) 



(30) 



and 



Ffn = 



e 2 \p 
Aire R 



F ±R (cf>) E(l-/?cos(<£/2))<# = 
/?(l-/3 2 ) 



(3 2 ln(tan(0/2)) 



l-/3cos(0/2) 



(3,1) 
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where the superscript "B" stands for "bunch", and Ao 
is the constant linear density of the bunch. If we now 
expand in </> m ; n <C 1 and ^ -C 1 the trigonometric 
functions in Eq. (^o|) and Eq. (|3l| ) we have 



Ft 



e 2 A 

47T£ -R 



— Ill 



2 7 2 V 1 - + /% ax /8 1-/3 + /30Ln/8 



(32) 



e 2 A 
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"2 ( As" - 



and 




AneoR 



? 

max 



- 2As 2 



(35) 



(36) 



and 



T 



F 



_LR 



e 2 A 
iireoR 



hi 



1 



1 - (3 + /3(f)'. 



max/ 



1-/3 + ^/8 



(33) 



The second term on the right side of Eq. (|32j) is centrifu- 
gal, as well as the first one on the right side of Eq. (33). 
The other ones are centripetal. Note that the logarith- 
mic term in Eq. (^) is unimportant (in the limit of large 
values for 7) with respect to the one in Eq. (|33|), while 
the second term in Eq. (32) modifies for a factor 1/2 the 
analogous centripetal term in Eq. (^). Therefore, the 
total transverse force on the test electron is given by: 




4ire R 



1 

~ T 

AS'; 



, 2 

max 



(37) 



From Eq. ([57|) we see that the centripetal term tends, 
asymptotically, to zero as (70 ma x) 2 - 

In the case </> max 3> I/7, instead, we have 



F 



_LC 



e 2 A 
At:e q R 



ATTEnR 




__L m ( ( 24A§ ™a* 

7 2 \ As 



(38) 



F 



e 2 A 
Ane R 




b 2 < 

min / 



(34) 

which is the sum of a logarithmic centrifugal term and a 
centripetal term. It is easy to check that Eq. ( |3^ ) can 
be also obtained by direct integration of Eq. (|l^) . 

A natural assumption is to consider our bunch density 
such that there are many particles within a distance R/j 3 
which, as it has already been underlined, is the minimal 
characteristic distance of the system: if the test particle 
is to be considered at the head of our bunch, then it 
is straightforward to assume As m i n <C R/y which, by 
means of the retardation condition Eq. (|16|), gives us 
back the condition: m i n <C 1/7, the non-linear term in 
ofEq. @ being, in this case, negligible. 

Under the latter assumption we can easily investigate 
the two cases of a short bunch, that is </> max < I/7 (in 
which the linear term in </> max of Eq. ( |l6| ) dominates), 
and of a long bunch, that is </> max 3> I/7 (in which the 
linear term in max of Eq. (|l6|) is negligible) . These two 
cases correspond, of course, to the asymptotic situations 
for the two-particle interaction discussed in Section O. 

Let us consider first the case max <C I/7. Eq. p2[), 
Eq. (H) and Eq. @, then, read 



F 



e 2 \p 

An£ R 





and 



4ne R 



F 



e 2 \p 
Atte R 



- 2 



- 2 



(39) 



e 2 A 
AtteqR 



In 



In 



/ (24As E 



a/3 



At 



(40) 



which means that the centripetal term saturates to a con- 
stant value in the limit of a long bunch. The case of a 
coasting beam (which fits our particular case ^ max ^ 1) 
with transverse extent h ^> R/j 2 has been discussed 
in literature (see p3[ , p[). As already mentioned 
above, the presence, in the choice of the model made in 
the latter works, of a non- negligible transverse dimension 
alters the structure of the retardation condition (with re- 
spect to the one given in Eq. (|l6|)), which explains why 
the centrifugal part of the transverse force scales, in liter- 
ature, as \n(R/h). Moreover, again because of the pres- 
ence of a transverse dimension, the velocity contribution 
to the transverse force can be discarded giving the same 
centripetal term which we found from the radiative term 
only in Eq. (|9}, that is -2e 2 A /(47r£ i?)- 
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Note that, in order to obtain Eq. (39), we had to inte- 
grate contributions from <f> <C 1 to <j> 1: we can give a 
simple explanation for the centripetal constant force term 
just analyzing Eq. ( |l3f) and its asymptotic behaviours. 
The product of Eq. (y|) with the Jacobian (Eq. (§|)) 
of the transformation between As and cj> (this product is 
just the integrand in Eq. |3l])) is equal, in the limits for 
<f> 3> 1 and <j> -C 1, to 1/0: then we can conclude that the 
logarithmic centrifugal term in Eq. (|39| ) (as well as in Eq. 
( |38| and Eq. (fto|)) takes into account the 1/0-behaviour 
of the transverse force for the two-particle system, while 
the constant centripetal term brings information about 
the way in which the transverse force for a two-particle 
system changes in the intermediate region between the 
limits <f> <C 1 and > 1. Note, however, that there is 
no physical ground to distinguish between the first (cen- 
trifugal) and the second (centripetal) term in Eq. (40): 



Let us designate with y the distance, along the straight 



the total force is always centrifugal, and both terms are 
consequences of the integration of a unique expression 
for the force between a two-particle system (which, of 
course, is always centrifugal too) . 



IV. TRANSVERSE INTERACTION BETWEEN 
TWO ELECTRONS MOVING IN AN ARC OF A 
CIRCLE 

We will now discuss the case of a two-particle system 
during the passage from a straight path to a circular one 
and from a circular path to a straight one. The four 
possible cases are sketched in Fig. |^ for the case in which 
the test particle is in front of the source. 

The case in which both particles are in the bend, de- 
picted in Fig. ^b, has already been discussed in Section 
O. Moreover, note that the situation in which the source 
particle is ahead of the test electron can be treated im- 
mediately for all three (a, c and d) transient cases in 
Fig. H (of course, with respect to the figure, test and 
source particle exchange roles) on the basis of Eq. (19). 
In fact in such a case, we can assume the retarded an- 
gle (j> small enough (the test particle "runs against" the 
electromagnetic signal) so that the actual trajectory fol- 
lowed by the particles is not essential and one can use Eq. 
( |l9| ) to describe also the transient cases. Now, the impor- 
tant contribution from the source particle comes from the 
acceleration part of the Lienard-Wiechert fields. Then, 
within our approximations, the only non-negligible con- 
tribution is constant and identical to the one in Eq. ( |l9| ) , 
and it is present in the situation (again, with the roles 
of test and source particle inverted) depicted in Fig. ||a 
only (Fig. |B|b is just the steady-state case). 

We will discuss more extensively the consequences of 
this fact in Section |v|. 

Let us now focus on the cases in which the source par- 
ticle is behind the test particle and, in particular, let us 
first deal with the case in Fig. |k; such a situation occurs 
when the following condition is met (see 0): 




FIG. 5: Relative configuration of the retarded source point 
S'and the test point T for a system of two electron passing a 
bending magnet. 



line before the bend, between the retarded source par- 
ticle and the beginning of the magnet. The retardation 
condition, in its approximated form, reads (see 0) 



Ay 



24 



(42) 



where we introduced the normalized quantity y = yj/R, 
which is just y/R normalized to the synchrotron radiation 
formation angle, I/7. 

In the situation considered, the source particle is only 
responsible for a velocity contribution, therefore F± = 
F±c- By direct use of Eq. (||), one can find the exact 
expression for F± 



As > 



24' 



(41) 
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Fi = 



(3 2 R{1 — cos 4>) — /3 sin cj)[(y + R sin ( 



R 2 (l 



COS( 



V^ + R-Rcost 



y sin < 



47T£o7 



{ [{y + R sin <f>) 2 + R 2 (1 - cos <p)} 1 / 2 - [3y - R/3 sin < 
I 



(43) 



Expanding the trigonometric functions in Eq. 
using normalized quantities one finds: 



3h and 



4 7 E 



+ <P? 



47re i? 2 
y 2 + y(^ + 3 /2) +0 4 /4 



(y + 0) 2 + ^ 4 /4 



(44) 



Eq. 



It can be easily verified that, as it must be, 
reduces to Eq. (^|) in the limit y — > 0. 

It is now possible, by means of Eq. (ff~4|), to plot the 
normalized transverse force F = i^/[e 2 /(47T£oAs)] as a 
function of the position after the injection (defined by the 
entrance of the test particle in the hard-edge magnet) for 
different values of As = Asj 3 / R. In Fig. ^| we com- 
pared such a plot with numerical results from the code 
TRAFIC 4 (see §). 

Note that, at the position which corresponds to the 
entrance of the retarded source in the magnet there is a 
discontinuity in the plots. This is linked to our model 
choice, and it is due to the abrupt (hard edge magnet) 
switching on of the acceleration fields. 

As general remark to Fig. || (and to the following ones) 
it might be worthwhile to stress that the perfect agree- 
ment (with graphical accuracy) between our calculations 
and numerical results by TRAFIC 4 provides, per se, an 
excellent cross-check between analysis and simulations, 
which enhance one's level of confidence on both these 
approaches. 



J 



Let us now consider the case depicted in Fig. |5|c, in 
which the source particle has its retarded position inside 
the bend and the test particle has its present position 
in the straight line following the magnet. We will de- 
fine with x the distance, along the straight line after the 
magnet, between the end of the bend and the present po- 
sition of the test particle. In this situation the following 
condition is verified (see j|]): 



As < 



+ x 



+ 4x 



24 



(45) 



where cp m = 70 m , (p m being the angular extension of 
magnet, and x = jx/R (the reason for this normalization 
choice for x is identical to that for y) . 
The retardation condition reads 



As 



ip + x ip 3 ip + 4.x 



(46) 



2 24 ^ + £ 
In contrast with the case of Fig. ||a, here we have contri- 
butions from both velocity and acceleration field. Again, 
by direct use of Eq. (|) and Eq. (|) one can find the 
exact expression for the transverse electromagnetic force 
performed by the source particle on the test particle 



F, =F, 



(47) 



where 



„2 R(l - cost/i)(1 - f3 2 cosip) - fismip ((x + Rsinip) 2 + R 2 (l - cosip) 2 ) 1/2 - (3x - pSRsmip 



47T£ 7 2 



and 



*±R 



Awen R 



1 /2 

((x + Rsinip) 2 + R 2 (l — cos ip) 2 ) 1 — fix cos ip — p3R sin ip 



■ ■ 1/2 

0i ( — ip ((x + R sin i/j) 2 + R 2 (I - cosip) 2 ) + fix + (3R sin ip - (3R(1 - cosip) s'mip 



n I \{{x + RsiniP) 2 + i? 2 (l - cosV>) 2 ) 1/2 ~ f3xcosip - /3Rsmip 

(R + x sin ip — R cos ip) 

1/2 

((x + Rsmtp) 2 + R 2 (l — costp) 2 ) ' — p3x cos ip — (3R sin ip 
R{\ -cosV0(l -0 1 cos ip*) -fismi); \((x + Rsin 1P) 2 + R 2 {\ - cosip) 2 ) 1 ' 2 - fJx - p3R sin ip 



(48) 



(49) 
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FIG. 6: Normalized transverse force (F = Fj_/[e 2 /(Atyso As)]) for a two-particle system entering a hard-edge bending magnet 
as a function of the position after injection. The solid lines show analytical results; the circles describe the outcome from 
TRAFIC 4 . We plotted several outcomes from different values of the normalized distance between the two particles. 



Expanding the trigonometric functions in Eq. (|4^) and Eq. (p^), and using normalized quantities one finds: 
I 



2e 2^3 -2.i 2 + x (V 3 - 2tp) + ?A 4 /: 

Fm= '—^(x + i)f i '- 



Aire R 2 



{x + xp) 2 + (^ 2 /4)(2.t + ^) : 



(50) 



and 



2eV 



x 2 + x 



V3(2 - t/. 2 ) + i, 2 - 3/4V> ( f + + ^/ 2 )^ h + ^(~ 2 + ^ ) + ^ V 2 



(i + V) 2 + (^ 2 /4)(2x + V)' 



(£ + V^) 2 + (^ 2 /4)(2i + ^) 2 



Simil arly to the latter case, it can be easily verified that respectively, in the limit x — > 0. Again, it is possible to 
Eq. ( pOl) and Eq. (^Tj) reduce to Eq. (g) and Eq. (10), plot the normalized transverse force F (defined above) as 
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Parameter: As y 3 /R 
1.0 2.0 5.0 




1.00 1.02 1.04 1.06 1.08 1.10 

Position after ejection (m) 
(ejection at s = 1 .0 m ) 

FIG. 7: Normalized transverse force (F — F±/[e 2 / (Aiveo As)]) for a two-particle system leaving a hard-edge bending magnet 
as a function of the position after the ejection. The solid lines show analytical results; the circles describe the outcome from 
TRAFIC 4 . We plotted several outcomes from different values of the normalized distance between the two particles. 



a function of the position after the ejection (denned by 
the exit of the test particle from the hard-edge magnet) 
for different values of As = As"/ 3 /R. In Fig. [7] we com- 
pared such a plot with numerical results from TRAFIC 4 . 

Again, at the position which corresponds to the exit of 
the retarded source from the magnet there is a discon- 
tinuity in the plots. This, again, is linked to our model 
choice, and it is due to the fact that, for particles on axis, 
when the retarded source leaves the magnet there is only 
Coulomb repulsion along the longitudinal direction. 

The last case left to discuss is depicted in Fig. ||d; the 
source particle has its retarded position in the straight 
line before the bend, and the test particle has its present 
position in the straight line following the magnet. This 
case occurs when 



As > 



+ 4i 



24 



Thc retardation condition reads 



As 



24 



(52) 



4>m + x + y 0f n <p m + 40 m (x + y) 



y 



(53) 



In this case we have only velocity contributions. The 
exact expression for the electromagnetic transverse force 
on the test particle is 



e 2 r r 
^\ (R + 2/sin<^ m — Rcos(j) m )(l — (3 2 cos^ m ) — (3 sin0 m (2R 2 + x 2 + y 2 + cos0 m (2a;y — 2R 2 ) 

47T£o7 I L v 
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+2R(x + y) sin m ) 1/2 - 0(x + ycos^ m + R sin (j> m ) } j [2R 2 + x 2 + y 2 + 2 cos m (xy - i? 2 ) + 2R(x + y)sm4> m ] 

1 -3 

— /? cos m (a; + y cos m + R sin </> m ) — /3 sin <f> m (R + y sin m — R cos ^> m ) > (54) 



,1/2 



r 



Expanding the trigonometric functions in Eq. (|5J) and proximated expression for F±: 
using normalized quantities one finds the following ap- 

I 



i 7 ! 



4ire a R 2 



8 7 3 (£ + ii + 0m) 2 ^ m {~x 2 /2 + f/2 + (4> 2 j2)xy + £(&/4 - 4> m /2) + y(4> 3 j4 + m /2) + ^f n /8} 

x{(x + y + 4> m ) [x(l + 4> 2 m )+y + 4> m + 4> 3 m /3\ - ($ 2 J12) [I2xy + 4(4 + y)4> m + &]} ' .(55) 



It is easy to verify that Eq. (p5|) reduces, respectively, 
to the steady state (Eq. (H^V) when x — and y = 0, 
to the transient case in Fig. ||a when x — (Eq. (0)) 
and to the transient case in Fig. |5|i when y = (Eq. 
(fl8|)). Following the treatment of the transient situations 
in Fig. |^a and in Fig. ||c it is possible to plot, for this 
case too, a normalized expression for the transient force, 
i.e. the usual F, as a function of the curvilinear position 
of the test particle (s=0 indicates the entrance of the 
magnet) for different values of As = As^/R and for 
different magnet lengths. In Fig. ||, Fig. |jj and Fig. [l(] 
we compared our analytical results with numerical results 
from TRAFIC 4 , for the cases As = 0.2, As = 1.0 and 
As = 5.0 respectively. 



V. TRANSVERSE INTERACTION BETWEEN 
AN ELECTRON AND A BUNCH ENTERING A 
BEND FROM A STRAIGHT PATH 

In the previous Section we dealt with all the possible 
configurations for a two-particle system moving in an arc 
of a circle. Now, one can consider a bunch moving on 
the same trajectory and calculate the transverse force on 
a test particle as the sum of contributions from all the 
source particles within the bunch. 

As an example, we will study here the case of a bunch 
entering a long bending magnet. Such a case is impor- 
tant, as mentioned before, for code benchmark purposes 
and for direct application in restricted regions of param- 
eters (negligible transverse bunch size and bunch energy 
spread) . 

First we will analyze the case of a bunch with rect- 



J 



r 



angular density function, and we will assume the test 
particle to be behind the bunch. Such an analysis will be 
performed using our line model and it is therefore valid 
only for transverse dimension of the bunch much smaller 
than the distance between the test particle and the bunch 
As m i n . After the discussion, in Section [V, about a two- 
particle system with the test particle behind the source 
electron, one is led to conclude that, within an electron 
bunch, interactions between sources in front of the test 
particle and the test particle itself are important and, in 
general, they must be responsible, at the entrance and at 
the exit of the bending magnet, for sharp changes in the 
transverse forces acting on the test electron. The quan- 
titative change depends, of course, on the position of the 
test particle inside the bunch: the extreme cases are for 
the test particle at the head of the bunch, where there 
are just interactions with electrons behind the test par- 
ticle (no head-tail interactions), and for the test particle 
at the tail of the bunch, where all the sources are in front 
of it (only head-tail interactions). It may be worthwhile 
to underline that the sharp jumps in the transverse force 
are expected to take place in a space interval compara- 
ble, at most, with half of the bunch length, in the case 
of the test electron at the tail of the bunch. In order to 
show this, one can easily calculate the transverse force 
acting on a test particle behind a bunch with rectangu- 
lar density distribution entering a hard-edge magnet. If, 
as usual, we indicate with As max the distance from the 
test particle to the head of the bunch and with As m i n the 
distance from the test particle to the tail of the bunch, 
then one can easily derive such an expression from Eq. 
<M) and Eq. ph: 



F± htWO - { e 2 A /(4TC i?) In [As„ 
I e 2 Ao/(47T£ i?) ln[As n 





x /(As E 
</As mi 



+ /?))] 



< 



■ < 
< (As„ 
> (As n 



■ As E 
As n 



a)/[A(l 

0/^(1- 



■0)] 
•/?)] 



(56) 
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FIG. 8: Normalized transverse force (F — F±/[e 2 /(47reoAs)]) for a two-particle system crossing a hard-edge bending magnet 
as a function of the position of the test particle inside the magnet in the case of a short magnet <j> m <C 1. The solid lines show 
analytical results; the circles describe the outcome from TRAFIC 4 . Here the normalized distance between the two particles is 
As = 0.2. We plotted several outcomes from different values of the magnet length R4>m = 0.1, 0.2, 0.3, 0.4, 0.5, 1.0. 



where "HT" stands for "head-tail". 

The following step is to actually plot the transverse 
force in Eq. (pq). It is convenient to choose, as nor- 
malization factor for the transverse force, the value f± = 
e 2 Ao/(4ir£oR) ln[As max /As m in]- Our results, compared, 
once again, with simulations by TRAFIC 4 , are shown in 
Fig. |ll] for a bunch length As max — As min = 100 /xm, 
7 = 100, R = 1 m and for As m i n = 1 fim. This means 
that, according to the validity limits of our model, this 
result can be applied only when /i<l /im. Nevertheless, 
this simple example shows that the code TRAFIC 4 is ac- 
tually able to account for head-tail interactions. Note 
that, as expected, the transient has a spatial extent of 
about one half of the bunch length (As max /(1 + (3)). 

We will now analyze the case of a bunch with rectan- 
gular density function with the test particle in front of 
the bunch, as depicted in Fig. |j. In the injection case 
we have contributions from retarded sources both in the 
bend and in the straight line before the bend. The contri- 



bution from the retarded sources in the magnet is given, 
basically, by Eq. (p4|), and reads 




where " m" reminds that the contributions treated by Eq. 
([57]) are all from the " magnet" . All that is left to do 
now, is the investigation of the values which (f) m - ln and 
^max assume. Let us first define with <p* the solution 
of the retardation equation As min = (j)*/2 + 0* 3 /24. If 
4>* < <j>, the retarded position of the first source particle 
is in the bending magnet, and </> m i n = <j>* ■ On the other 
hand, when <ff > <j> there are no contributions to the 
transverse force from the bend. Next, we define with 
</>** the solution of As max = </>** /2 + 0** 3 /24. Supposing 
4>* < (f>, if </>** < 4> too, then all the particles contribute 
from the bend, and c6 max = df* . On the other hand, when 
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FIG. 9: Normalized transverse force (F — F±/[e 2 /(47reoAs)]) for a two-particle system crossing a hard-edge bending magnet 
as a function of the position of the test particle inside the magnet in the case of a short magnet <j> m <C 1. The solid lines show 
analytical results; the circles describe the outcome from TRAFIC 4 . Here the normalized distance between the two particles is 
As = 1.0. We plotted several outcomes from different values of the magnet length R4>m = 0.1, 0.2, 0.5, 1.0, 1.5, 2.0. 



4>** > <[>, we have a mixed situation, in which part of the 
particles contribute from the bend and others from the 
straight line before the magnet. In this case max = 4>- 

The contribution from the retarded sources in the 
straight path before the bend is given by 



As„ 



*F x (y(As),4>)dAs, (58) 



stands for " straight path" , and where the ex- 
pression for F± in the integrand is given by Eq. (44). It 
is convenient, as done before, to switch the integration 
variable from As to y. The Jacobian of the transforma- 
tion is given by (see 0) 



dAs 
dy 



2(0 + y) 2 



(59) 



After substitution of Eq. (|5J) and Eq. (|4J) in Eq. (|58|) 
one can easily carry out the integration, thus getting 



ft 



2e 2 A 



+ 20 + <fr 



4y min + 8y min + 40 2 +0 4 

(4j/ max + 2<p + (f> 3 



4j/i 



,",2 
max 



t </> + W + <p A 



(60) 



As done before for </> m i n and m ax, we can now investigate 
the values of y m in and y max . Let us start with y m in- First, 
we define with y* the solution of the retardation condi- 
tion AJ min = + y*)/2 + (0* 3 /24)(4y* + + $). 

If y* > 0, the retarded position of the first source par- 
ticle is in the straight line before bending magnet, and 
ymin = V* ■ On the other hand, when y* < 0, the re- 
tarded position of the first source particle is in the bend, 
and y min = 0. 

Next, we define with y** the solution of As max = 
(4> + y*)/2 + (0* 3 /24)(4y* + <j>)/(y* + </>). Consider the 
case y** < 0: all the particles contribute from the bend, 
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FIG. 10: Normalized transverse force (F = _Fx/[e 2 /(47reoAs)]) for a two-particle system crossing a hard-edge bending magnet 
as a function of the position of the test particle inside the magnet in the case of a short magnet m <C l.The solid lines show 
analytical results; the circles describe the outcome from TRAFIC 4 . Here the normalized distance between the two particles is 
As = 5.0. We plotted several outcomes from different values of the magnet length Rcj) m = 0.5, 1.0, 2.0, 3.0, 4.0, 5.0. 



that is we entered the steady-state situation. In this case 
2/max = ?/min = 0. On the other hand, when y** > 0, 
we have again a mixed situation, in which part of the 
particles contribute from the bend and others from the 
straight line before the magnet. In this case y max = y** ■ 

The following step is to actually plot the transverse 
force on an electron from a bunch with rectangular dis- 
tribution entering a long bend. It is convenient to choose, 
as normalization factor for the transverse force, the value 
/ = e 2 Ao/(47T£o-R) ln(As max ). Our results, compared, 
once again, with simulations by TRAFIC 4 , are shown in 
Fig. |l| for a bunch length of 100 fim, 7 = 100, R = 1 m 
and for different values of As. 

It is interesting to calculate the asymptotic expression 
for F± in the limit for a long bunch (As max 3> 1) and for 
a short distance between the test particle and the head of 
the bunch (As m j n «C 1). First let us indicate with 0b the 
normalized angular extension of the bunch (0b = 70b, 
where 0b is just the ratio between the bunch length and 
the radius of the circle R) . By means Eq. it is easy 



to prove that, when 0/0b < 2~ 2 / 3 , y max takes bigger and 
bigger values, with an upper limit y max = 2As max . On 
the other hand, when 0/0b > 2~ 2 / 3 , y max takes always 
decreasing values, with a lower limit y max = 0. Assuming 
2/min *C 1, one may check that, in the long bunch limit, 
Eq. <M) reads 



e 2 A /(27re i?) 




> b < 2- 2 / 3 
> b > 2- 2 / 3 



(61) 



which is a boxcar function. Note that in the passage from 
Eq. @ to Eq. @ we used the fact that <p > 1. In 
order to visualize the limiting process we plotted, in Fig. 
O, F® a , as it is given in Eq. ( |60| ) and normalized to / = 
e 2 A /(47re i?) ln(As max ) (i.e. v = F^Jf, in the plot), as 
a function of the position after injection, normalized to 0b 
(i.e. u = i?0/(70b), in the plot). The plots in Fig. 13 1, 
b, c and d refer to different bunch lengths (respectively 
10 2 /im, 10 3 ^m, 10 4 /zm and 10 5 ^m). For every choice of 
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FIG. 11: Normalized transverse force F±/fi acting on a test particle behind a bunch with rectangular density distribution. As 
the bunch enters a hard-edge bending magnet we plot the normalized force as a function of the position of the test particle inside 
the magnet (solid line) and we compare with results by TRAFIC 4 (circles). Parameters choice are: bunch length = 100/im, 
7 = 100, R = lm; the value of As m i n is 1/im. 



the bunch length we show different outcomes for several 
choices of As m ; n . As one can see from Fig. [l3], in the 
limit for As m i n -C 1 and As max 3> 1, one approaches 
the boxcar function described by Eq. (|6l"|). Note that, 
here, the radius of the bend is R = lm and 7 = 100, 
therefore we have 2 _2 / 3 7?/7 ~ 6.3 • 10~ 3 m, while the 
maximum value for v, due to the normalization choice, 
is given, from Eq. (pl|) , by 2/ln(As m i n ). In the case of 
Fig. |l3|d, for example, we have 2/ln(As m i n ) ~ 0.174 (in 
agreement, of course, with the maximum value found in 
the plot). 

On the other hand, again in the same limits, one can 
find the asymptotic expression for the contribution from 
particles with retarded position in the bend: this is given, 
indeed, by Eq. (|57|), with the limit expressions for </> min 
in the case of a small distance As m j n , and max in the 
long bunch limit. For <^ m ; n we have, from the retardation 
condition in Eq. (42) <fi m - nl = 2s m ; n . 



For </> max , if As max < 3 /24 we have, again from Eq. 
©, max = 0b = (24As max ) 1 / 3 . Otherwise, when 



As max > c/) 3 /24 then (f) max — <f>, and we have the steady 
state situation. 

By means of these results one can build an expression 
for the transverse force in the case of a bunch with general 
density distribution A, by considering it as a composition 
of rectangular bunches with length (s — s') and linear 
density ds'dX(s')/ds' , under the constraint that the im- 
portant bunches for such a composition are long enough 
to neglect all linear terms in the expression for the retar- 
dation condition (i.e. X(s)j 3 /R 3> d\(s')/ds', see gj). In 
this case, the contribution from particles in the straight 
line reads 



r s-R4> /6 



-ds 



-X(s - i?0 3 /6) (62) 



2ire R 

Note that the latter expression does not depend on 
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FIG. 12: Normalized transverse force (F±/f) acting on a test particle from a bunch with rectangular density distribution 
entering a hard-edge bending magnet as a function of the position of the test particle inside the magnet. The solid lines show 
analytical results; the circles describe the outcome from TRAFIC 4 . We chose As max = 100/im, 7 = 100, JJ = lm; graphs are 
plotted for several values of the parameter As m i n . 



As m in- On the other hand the contribution from par- 
ticles in the bend is 



j^tot 
^_Lm 



S - s') 



,.d\(s') , 



ds 1 



ds 



27T£ i? J- 



1, /24 7 3 A 

- 1 + 3 ln hr AsE 



— ds , (63) 



which, instead, depends on As m j n . It is very interesting 
to show that this dependence on As m i n cancels with the 
dependence on As m i n of the steady-state force: in fact 
this constitutes a general result independent from the 
choice of the position of the test particle. In order to 
show this, let us first note that Eq. (p3|) can be written 



as 



F 



±m 




1 +ln(<£b(s')) -ln(^ min ) 



-1 + In(^) - ln^min) 



4:TTS R 

ds 
ds' 



(64) 



while the steady state contribution is given by 

e 2 

ptot ^ c x 

Xstcady - 4ir£oR 

-1 + H4> h (s')) - ln(0 min )] ^r-ds' .(65) 

Subtracting side by side Eq. ( |65|) from Eq. and 
adding the contribution from the straight path, one fi- 
nally gets the following "regularized" expression for the 
transient transverse force: 
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FIG. 13: Plot of v = Ff s // as a function of u = R^/^cj)^); here i? = lm and 7 = 100. Case (a) the bunch length is 10 /im, 
(b) 10 3 pm, (c) 10 4 /im and (d) 10 5 /im. Each case is parameterized with respect to different values of As m i n = 7 3 As m i n /_R (see 
the legends in the plots). 



VI. TRANSVERSE DYNAMICS AND 
CANCELLATION OF CENTRIFUGAL FORCES 



jptot 



s-R<j> 3 /24 



hi 



2we R 
24(a- 



X(s - R(j> 3 /6) 
s')\ d\{s') 



ds> 



ds' 



(66) 



which is completely independent from As m ; n and, there- 
fore, free from singularity in the limit As m ; n — > 0. It 
might be worthwhile to remark that usual regularization 
techniques take place, in the study of longitudinal (CSR) 
self- interactions (see [||) at the stage of the two-particle 
system and before integration of the contributions from 
all the retarded sources within the bunch. The situation 
is reversed here, where regularization takes place after 
integration. 



We will discuss here, from a qualitative viewpoint, the 
effect of perpendicular forces on the transverse beam dy- 
namics (although, as we stated in Section Q, this article is 
mainly devoted to the study of electrodynamical effects). 
In particular we will concentrate on the controversial is- 
sue of the cancellation of the centrifugal force, in the 
transverse equation of motion of a particle, by the ratio 
between the electron energy deviation of from its nom- 
inal value and the design radius R (sec M, H, 

H). 

To illustrate this cancellation as it is explained in liter- 
ature, let us consider a bunch moving in a circular orbit 
with design radius R, and let us indicate with X the 
transverse displacement of a test electron from the equi- 
librium orbit. Then (see, for example, H), up to the 
first order in X , one can easily write the equation for the 
transverse motion of the test electron as: 



20 



d 2 X 
c 2 dt 2 



X 
R 2 



AEo 
RE 



(67) 



where Eq — 70 mc 2 is the design energy (which is linked 
to the equilibrium radius by the relation /3~/amc = eBR), 
AEq = (7 — 7o)toc 2 is the kinetic energy deviation, from 
Eq and, finally, F r is the self-interaction force in the 
transverse direction. 



We already discussed (see Section III) the fact that the 
transverse force on a test electron, both on-orbit (X = 0) 
and off-orbit (X ^ 0) can be written, in the long bunch 
limit, as the sum of two terms: a logarithmic, centrifugal 
term, and a constant, centripetal one. 

In jl8| , the case of bunch compression by means of a 
magnetic chicane (at LCLS) is analyzed, and the cen- 
trifugal term is told to be essentially (aside for a neg- 
ligible residual) canceled by part of the Ai?o-term in 
Eq. (|67]). In fact all contributions to the transverse 
emittance growth are attributed to the centripetal force 
" which originates from radiation of trailing particles and 
depends on the local charge density along the bunch. The 
maximum force takes place at the center of the bunch and 
its effect on the transverse emittance is estimated in the 
reference. This estimate predicts an emittance growth of 
■C 1% for the worst case (last dipole of chicane-2 where 
the bunch is shortest)" (quoted from |fl8|). 

Here we will discuss this effect, well known to the CSR 
community as "the cancellation effect", within the limits 
of our line model, and show that, in contrast to what has 
been implicitly assumed in [ fl8| , it has no general validity. 

As a general remark to all the previous analysis of the 
problem, we must say that, up to now, only the situation 
of sources behind the test particle has been discussed, 
while we know, from Section IV, that head-tail interac- 



tions are present too: they are characterized, at least 
within the limits of our model, by a magnitude of the 
same order of the tail-head interactions and the spatial 
extension of their transient is negligible (comparable, at 
most, with half of the bunch length). Since this centrifu- 
gal force depends on the position of the test particle along 
the bunch, it will be responsible for normalized emittance 
growth. 

Anyway, let us consider, in particular, the tail-head 
interaction problem. As underlined in Section |, the full 
solution to the evolution problem is met when one is able 
to solve simultaneously the equation of motion and the 
equation for the electromagnetic field. Then, in princi- 
ple, we may adopt, in our discussion, two separate view- 
points: in the first one we imagine to solve, in some way, 
the full evolution problem, while in the second we treat 
the same situation by means of a perturbation theory 
approach, assuming that, in first approximation, the mo- 
tion of a rigid bunch is driven by the external magnetic 
field alone and then calculating the perturbation to the 
particle motion due to transverse self- fields: we will show 
that, in both cases, the cancellation effect has no general 
validity. 



Let us begin with the first viewpoint. Generally speak- 
ing, if there was complete compensation between the to- 
tal transverse force and the A-Eo-term, then the particles 
with the same total energy (assumed equal to the sum of 
kinetic energy and potential self-energy from the bunch) 
would have followed the same trajectory. Such a com- 
plete cancellation would be based upon two assumptions. 
Firstly, AE = -eA$ and, secondly, F± = e$/i?. Note 
that, in this case, AE + RF± = e$(0). This, of course, 
is not the case, since the cancellation has always been 
understood for the centrifugal term alone (and, anyway, 
it is not complete). Moreover let us remind that, as ex- 
plained in Section |TTl| , there is no physical ground to dis- 
tinguish between the centripetal term and the centrifugal 
one: from a physical viewpoint there is just an overall 
centrifugal force. A further subdivision is just of math- 
ematical nature, and explains how the transverse force 
plotted in Fig. ^ behaves from one asymptotic (short 
distance between test and source particle) to the other 
(large distance between test and source particle). This 
fact, alone, suggests that the cancellation issue has no 
general validity and that, indeed, it is artificial (as the 
subdivision between centripetal and centrifugal term is) 
even when it works, like in the case of a coasting beam 
in a simple, circular steady-state motion. 

Nevertheless, let us retain, in our terminology, the dis- 
tinction between centripetal and centrifugal term. Then, 
a complete cancellation between centrifugal force and po- 
tential term would be based, similarly as before, upon two 
facts. Firstly, AE = -eA$ and, secondly, i^cntrifugai = 
e$/i?, at any time. Whenever one of these assumptions 
is no more verified, the cancellation fails to happen. 



Of course, in order to show that the cancellation ef- 
fect is not valid in general, it is sufficient to provide a 
counter-example. Let us consider, therefore, the case of 
a bunch with zero initial energy spread or, more simply, a 
two-particle system with the test electron in front of the 
source, such that A7 = 7 — 70 = 0. Let us restrict to a 
case in which only space charge effects are important, as 
concerns the longitudinal motion: therefore we consider, 
in this sense, two electrons running on a straight path. 

The Coulomb (space-charge) interaction changes both 
the kinetic energy of the particles and their relative dis- 
tance, which we will assume, initially, equal to As in 
the laboratory frame. Let us deal with the problem in 
the (instantaneous) rest frame of the center of mass of 
the system (simply designated as the rest frame, in the 
following) and assume that the particles do not move 
relativistically in such a frame. After a long time one 
will be left, asymptotically, with two particles far away 
from each other, each one with a kinetic energy equal to 
e 2 /(87reoAs), where As = 7AS. Therefore, the change 
in kinetic energy of the test particle in the laboratory 
frame, can be found by means of a Lorentz transforma- 
tion (performed on the momentum in the longitudinal 
direction) : 
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. c ( 2me 2 \ c ( ^me 2 \ ^ 2 , . 

AEo =^{s^As) = p{^As) • (68) 

On the other hand, again in the asymptotic limit, the 
potential energy of the test electron (as well as the one 
of the source) will approach zero. Then, at least after a 
long time from the beginning of the evolution, one has 

AE (t -> oo) ^ -eA$(i -» oo) = e$(0) . (69) 

This means that the first assumption for the cancellation 
(AE Q = — eA$ at every time) is violated, due to the in- 
clusion of space charge forces in the longitudinal motion; 
as a result we can say that there is no cancellation effect. 

Of course, one could also treat the situation of a circu- 
lar motion, which is more complicated, but the example 
we gave is sufficient to show that the effect is not valid 
in general. 

We are now left to deal with the second viewpoint. 
As it has already been said above, by following this ap- 
proach, we assume that, in zero approximation, the par- 
ticles move under the action of the external field alone 
and then, once the particles motion is fixed (in the ze- 
roth order), one can calculate the perturbation to such a 
motion due to self-fields. Note that, if one treats in this 
framework the example of a coasting beam or of a rigid 
line bunch with finite longitudinal dimension, then $ is 
constant, which suggest there are no more contradictions 
with the cancellation assumptions. Nevertheless we will 
now show that cancellation is far from being a general 
effect. Let us study, first, the case of a rigid bunch. Con- 
sider, for example, the two cases of a rigid line bunch 
oriented either along the velocity or perpendicular to it 
(see Fig. m . 

For simplicity we can study the non-relativistic case: 
in the relativistic one we have an analogous situation. 
The kinetic energy of the line bunch and the scalar po- 
tential acting on each particle do not depend on the ori- 
entation of the bunch with respect to the velocity vector. 
If the cancellation effect had general validity, also the 
transverse self-force should not depend on the orienta- 
tion but, indeed, one can easily verify that, for the case 
of the line bunch oriented along the perpendicular direc- 
tion (see Fig. |l4|b) , the self- interaction force is two times 
smaller than in the case of orientation along the longitu- 
dinal direction. In fact, in the case of the bunch oriented 
in the longitudinal direction we have repelling forces par- 
allel to the velocity vector while, in the other case, they 
are perpendicular: according to Eq. (^3|), the momen- 
tum flux along the parallel direction (that is T33) is then 
equal to the potential energy but having opposite sign in 
the first case, and equal to zero in the second. This fact 
shows that the cancellation effect has no general validity. 
This can be discussed in a somewhat deeper way too. 
As we showed in Section H, in order to get a correct dy- 
namical interpretation of the electrodynamical transverse 




O (b) 



FIG. 14: Schematic of a line bunch in a circular trajectory; 
(a) the line bunch is oriented along the velocity; (b) the line 
bunch is oriented perpendicularly with respect to the velocity. 



forces for a two-particle system, one needs to account for 
the self-interaction energy and momentum flux and, by 
doing so, the energy and momentum of the system do 
not form a four-vector anymore (as, instead, people al- 
ways assume in considering the cancellation issue) . They 
are, indeed, components of a second rank tensor (the 
energy-momentum tensor) and, therefore, transform in a 
different way with respect to a four- vector: this is the 
conceptual reason which explains why the cancellation is 
not a fundamental effect; as soon as one cannot consider 
energy and momentum as a four- vector anymore the can- 
cellation is spoiled, in general, by the T33 component of 
the stress tensor in the rest frame. 

Also the analysis of the transient behavior leads to the 
conclusion that the compensation is not of general im- 
portance. Let us consider the case of a bunch entering 
a bend (which we discussed in Section M, in the labo- 
ratory frame. The formation length for both the scalar 
potential and the transverse force (after which the bunch 
reaches the steady state situation) is simply the over- 
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taking length 2(3si? 2 ) 1 / 3 . Now, on the one hand the 
transverse force is zero at the beginning of the bend, and 
relaxes to a 7-independent value in the steady-state situa- 
tion while, on the other hand, the scalar potential starts 
from a 7-dependent value before the bend and relaxes 
to a 7-independent value in the steady state regime too. 
Comparing the two transients for the scalar potential and 
for the transverse force, we can therefore conclude that 
the first one is characterized by an additional free param- 
eter (the initial value of the beam energy) with respect 
to the second, and this means that the two transients are 
completely independent. As a result, again, the compen- 
sation proves to fail in this case (since it should follow 
from the subtraction of the centrifugal force with the 
term, in AEq/(REo), proportional to e<E>/i?). Of course, 
again, from a more general viewpoint, the failure of the 
cancellation can be just seen as a consequence of the fact 
that energy and momentum are components of a tensor, 
and not of a vector, that is as a consequence of the dif- 
ferent geometrical nature of a vector with respect to a 
tensor. 

As a last example, one can study, again in the per- 
turbative framework, the case of a beam with an initial 
kinetic energy chirp (which encompasses the bunch com- 
pression case) or, equivalently, the case of two electrons 
moving on rectilinear trajectory in such a way that that, 
for the test particle, 



AT, 



chirp 



A7(0)mc 2 > 



7me 
(3 V47re As 



1/2 



(70) 



which justifies the perturbative approach. 

It is possible to demonstrate that, in first order in the 
self-fields, space charge induces a change in kinetic en- 
ergy, in the laboratory frame, given by 



AEn = 



777ie 



1/2 



/AT, 



-(- , 

(3 V47re As J \ 87re As 



chirps 



1/2 



/AT c hirp 



where we started from the rest frame, we used a first or- 
der expansion in the kinetic energy difference, AT cri i r p, 
and, finally, a Lorentz transformation to the laboratory 
frame. Eq. ( f7l| ) shows that AE includes a free parame- 
ter (equal to (AT c hi rp e 2 /87reoAs) 1 ^ 2 /AT c hirp) which de- 
pends on the energy chirp. This makes the assumption 
ATo = — eA$ incorrect and demonstrates that the can- 
cellation issue does not exist at all during the bunching 
process. 

To sum up, the final result of this investigation is that 
the cancellation effect has by no means general validity 
and that one must exercise extreme care when dealing 
with transverse dynamics issues. For example, the qual- 
itative discussion in this Section shows that, in contrast 
with 18 , when one considers the problem of calculating 



the emittance dilution in a bunch compressor, not only 
centripetal force but also centrifugal contributions must 



be taken into account. We can conclude that, in prac- 
tice, the total force one has to account for is, at least, one 
order of magnitude larger than the so called centripetal 
force. 



VII. SUMMARY AND CONCLUSIONS 

In this paper we presented a fully electrodynamical 
study of transverse self-forces within an electron bunch 
moving in an arc of a circle. Our analysis is based on a 
line-bunch model. In the case of test particles in front of 
the source, our model is valid whenever h <C -R/7 2 . On 
the other hand, when the source electron is in front of the 
test electron, then the situation is more complicated and 
the model can be applied for |A|s ^ h The cases that do 
not fulfill these conditions are left for future study. 

In Section [n] we studied, first, the situation of a two- 
particle system moving on a circular path, finding an ap- 
proximated expression for the transverse self-forces which 
is the product of a factor dependent on the parameters 
that specify the system setup and a universal function 
(i.e. a function independent of such parameters): for 
every distance between the particles we found a centrifu- 
gal force, which has a qualitative and quantitative ex- 
planation by means of simple arguments from relativistic 
dynamics. We concluded that both the tail-head and 
head-tail interactions are important: in the first case the 
Coulomb interaction plays a role besides the Radiative 
one, while in the second only Radiative contributions are 
present. 



Further on, in Section [II , after discussing the applica- 
bility region of our model, we integrated the results for 
a two-particle system, thus finding an expression for the 
transverse interaction between a line bunch and a parti- 
cle in front of it. Such an expression is structured as the 
sum of a centrifugal logarithmic part and a centripetal 
term, of which we studied the asymptotic behaviors in 
the limit of short and long bunch, both with a small dis- 
tance between the test particle and the bunch head. 

In particular we found that, in the limit of a short 
bunch, the centripetal force tends to zero as (7</>max) 2 , 
while in the limit of a long bunch we found results already 
well known in literature. We were able to explain the 
constant centripetal term in the latter case as an overall 
effect of the transient between the asymptotic behavior 
(identical in 0, different in As) of the transverse force 
in the two-particle system, respectively for small or long 
distance between the two electrons. We concluded that 
the centripetal term is therefore only of mathematical 
nature, and there is no physical ground to distinguish 
it from the centrifugal term: from a physical viewpoint 
there is, in fact, just an overall centrifugal force. 



In Section IV, again within the region of applicability 



of our model, we extended our considerations to a two- 
particle system moving in an arc of a circle, thus finding, 
for the first time, exact and approximated analytical ex- 
pressions for the transverse force in all the possible tran- 
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sient configurations (see Fig. ||), including the case in 
which the source particle is in front of the test particle. 
Furthermore we plotted the expression for the transverse 
force in several practical cases, which are important for 
a quick evaluation of the magnitude of the effect and 
for cross-checks between computer codes. In particular 
we report a very good agreement with TRAFIC 4 , which 
demonstrates that such a code can easily deal, from a nu- 
merical viewpoint with the transverse transient problem. 

In Section [v|, we analyzed the situation of the trans- 
verse interaction between a line bunch and a test particle 
moving in an arc of a circle. We treated, in particular 
the case of the injection from a straight section into a 
hard-edge bending magnet. Firstly we calculated exact 
and approximated expressions for the transverse force; 
secondly, following what we did in Section [V, we pro- 



vided a few graphical examples; thirdly we analyzed our 
expressions in the limit for a long bunch and a short min- 
imal distance between the test particle and the head of 
the bunch. 

We showed that the contribution from the particles 
whose retarded positions are in the straight line before 
the bend is well described, as a function of the normal- 
ized angular position of the test particle inside the bend, 
by a boxcar function with length equal to the normal- 
ized bunch angular length multiplied by a characteristic 
constant which depends on the structure of the retar- 
dation condition. By simple composition of rectangular 
bunches we provided an expression for the calculation 
of the transverse interaction in the case of a long bunch 
with an arbitrary density distribution. We showed that, 



in the chosen limits, the contribution from the particles 
in the bend is independent from As m j n and that, in con- 
trast to this, the contribution from the particles in the 
straight line is dependent on As m ; n . Finally we proved 
that such a dependence can be removed by subtraction 
of the steady-state transverse self-interaction, thus pro- 
viding a "regularized" expression for F±. 

The case of injection provides a useful example for a 
quick evaluation of the transverse self-fields magnitude as 
well as for computer codes benchmark. By means of the 
same approach one can analyze also the case of ejection, 
which is left to future work. 



Finally, in Section VI we proved that the partial com- 
pensation of transverse self-force has by no means gen- 
eral validity, and one must exercise extreme care when 
dealing with transverse dynamic issues. In fact, in con- 
trast with fl8[| , we concluded that, when one considers the 
problem of calculating the emittance dilution in a bunch 
compressor, not only centripetal force but also centrifu- 
gal contributions must be taken into account. Therefore, 
in practice, the total force one has to account for is, at 
least, one order of magnitude larger than the so called 
centripetal force. 
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